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We obtained the analytic solutions of relic gravitational waves (RGWs) for the slow-roll inflation 
with a power-law form potential of the scalar field, V = \<f) n . Based on a reasonable range of n 
constrained by cosmic microwave background (CMB) observations, we give tight constraints of the 
tensor-to-scalar ratio r and the inflation expansion index fi for the fixed scalar spectral index n s . 
Even though, the spectrum of RGWs in low frequencies is hardly depends on any parameters, the 
»^ , high frequency parts will be affected by several parameters, such as n s , the reheating temperature 

■ Trh and the index fi s describing the expansion from the end of inflation to the reheating process. 
We analyzed in detail all the factors which would affect the spectrum of RGWs in high frequencies 

| including the quantum normalization. We found that the future GW detectors SKA, eLISA, BBO 

r ^ . and DECIGO are promising to catch the signals of RGWs. Furthermore, BBO and DECIGO have 

l"H ■ the potential not only to distinguish the spectra with different parameters but also to examine the 

£SJ ' validity of the quantum normalization. 
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W)| I. INTRODUCTION 

The validity of general relativity and quantum mechanics make sure the generation of a stochastic 
,_i ! background of relic gravitational waves (RGWs) [1-6] during the early inflationary stage, whose 
\ primordial amplitude could be determined by the quantum normalization at the time of the wave 

■ modes crossing the horizon during the inflation. Since the interaction of RGWs with other cosmic 

^ ; . components is very weak, the evolution of RGWs are mainly determined by the behavior of cosmic 

Q ! expansion including the current acceleration [7, 8]. Therefore, RGWs could serve as an unique tool 
^ ; to study the very early Universe earlier than the recombination stage when the cosmic microwave 
T"! ■ background (CMB) radiation was generated. As an interesting source for gravitational wave (GW) 
. ^ ■ detectors, RGWs exist everywhere and anytime unlike GWs radiated by usual astrophysical process. 
^ \ Moreover, RGWs spread a very broad range of frequency, 10™ 18 — 10 10 Hz, making themselves become 

one of the major scientific goals of various GW detectors with different response frequency bands. 
The current and planed GW detectors contain the ground-based interferometers, such as LIGO [9], 
Advanced LIGO [10, 11], VIRGO [12, 13], GEO [14], KAGRA [15] and ET [16, 17] aiming at the 
frequency range 10 2 — 10 3 Hz; the space interferometers, such as the future eLISA/NGO [18, 19] which 
is sensitive in the frequency range 10" 4 - KT 1 Hz, BBO [20, 21] and DECIGO [22, 23] which both 
are sensitive in the frequency range 0.1 — 10 Hz; and the pulsar timing array, such as PPTA [24, 25] 
and the planned SKA [26] with the frequency window 10~ 9 — 10 -6 Hz. Besides, there some potential 
very-high-frequency GW detectors, such as the waveguide detector [27], the proposed gaussian maser 
beam detector around GHz [28], and the 100 MHz detector with a pair of 75-cm baseline synchronous 
recycling interferometers [29] . Furthermore, the very low frequency portion of RGWs also contribute 
to the anisotropies and polarizations of CMB [30], yielding a magnetic type polarization of CMB as 
a distinguished signal of RGWs. WMAP [31-34], Planck [35], the ground-based ACTPol [36] and 
the proposed CMBpol [37] are of this type. 
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The reheating temperature, Trh, carries rich information of the early Universe, and relates to the 
decay rate of the inflation as T RH oc VT [38, 39]. Recently, the temperature of the reheating, T RH , 
was evaluated [40] according to the CMB observations by WMAP 7 [34], combining with the slow-roll 
inflation scenarios. Furthermore, the resultant RGWs was studied in [41]. However, these pieces of 
work underwent the assumption of a fixed form of the potential of the scalar filed driven the slow-roll 
inflation, V((p) = |m 2 2 . In this paper, we study a more general case of V((f>) = \cj) n [42], where A is 
constant. Moreover, we adopt the limitation n < 2.1 obtained from the spectrum of CMB [43]. For 
a non-fixed value of n, it is hard to evaluate the temperature of the reheating process, Trh? using 
the method employed in [40]. Thus, we choose several values of Trh lying in the range of ~ 10 4 — 10 8 
GeV, where the lower limit and the upper limit of Trh are obtained from the constraints in [43] and 
[44], respectively. In the text, one will see that n and Trh affect the increases of the scale factor in 
the early stages of the Universe. Once all the expansion histories of different stages are determined, 
the evolutions of the RGWs at various phases can be determined subsequently. For the present time, 
the solutions of RGWs can be obtained, whose different frequency bands correspond to the k-modes 
re-entered the horizon at different phases. On the other hand, the anisotropies due to the tensor 
metric perturbations (gravitational waves) can be scaled to those due to the observations of the 
scalar perturbations by introducing a parameter r called tensor-to-scalar ratio. Under the frame of 
the slow-roll inflation scenario, r will be constrained in a narrow range due to the constraints from n 
for a given value of the scalar spectral index n s . Similarly, the inflation expansion index will also 
be constrained in a narrow range. Besides, there is a simple relation between n and the preheating 
expansion index /3 S describing the expansion behavior of the universe from the end of inflation to the 
reheating process. As will be shown below, the RGWs in the very high frequencies are sensitively 
dependent on the parameters n s , j3 8 and T RH . Furthermore, the spectra of RGWs also depends on 
the condition of the quantum normalization. To this end, the spectra of RGWs given by different 
parameters and different conditions will confront the various current and planed GW detectors. The 
future detectors BBO and DECIGO are quite promising not only to determine various parameters 
but also to examine the validity of the quantum normalization. 

Throughout this paper, we use the units c = h = k B = 1. Indices A, /i, v,... run from to 3, and 
i, j, k,... run from 1 to 3. 

II. RGWS IN THE ACCELERATING UNIVERSE 

In a spatially flat universe, the existence of perturbations modifies the Friedmann-Robertson- 
Walker metric to be 

ds 2 = a 2 (r)[-dr 2 + (5ij + hi^cbfdxf], (1) 

where a(r) is the scale factor, r is the conformal time, and hij stands for the perturbations to the 
homogenous and isotropic spacetime background. In general, there are three kinds of perturbations: 
scalar perturbation, vectorial perturbation and tensorial perturbation. In this paper we only consider 
the tensorial perturbation, that is, gravitational waves. In the transverse-traceless (TT) gauge, 
satisfies: -g^f = and h l i = 0, where we used the Einstein summation convention. In the Fourier 
fc-modes space, it can be written as 

where a = +, x stands for the two polarization states, the comoving wave number k is related with 
the wave vector k by k = (tfyW) 1 / 2 , h£l*(r) = h[ a) (r) ensuring be real, and the polarization 
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tensor e\j satisfies [2]: 

$ ) e™ = 26 <fa ,, 4^ = 0, e£W = 0, e^(-k) = 4 CT) (k). (3) 
In terms of the mode the wave equation is 

h^'(r) + 2^h^(r) + k^\r)=0, (4) 

where a prime means taking derivative with respect to r. The two polarizations of h^\r) have the 
same statistical properties and give equal contributions to the unpolarized RGWs background, so 
the super index (er) can be dropped. The approximate solutions of Eq. (4) are well analyzed in 
[2, 3, 7], and are detailed listed in [41] given an accelerating universe at present. Furthermore, the 
analytic solutions were also studied by many authors [8, 45-48]. For a power-law form of the scale 
factor a(r) oc r a , the analytic solution to Eq.(4) is a linear combination of Bessel and Neumann 
functions 

h(r) = T^faJ^W + CtN^iikr)], (5) 

where the constants C\ and C 2 for each stage are determined by the continuities of h^r) and h' k (r) at 
the joining points ti,t s , r 2 and te [8, 45]. Therefore, the all the constants in the solutions of RGWs 
can be completely fixed, once the initial condition is given. In a spatially flat (k = 0) universe, the 
scale factor indeed has a power-law form in various stages [2, 41, 45, 47]. It is described by the 
following successive stages : 
The inflationary stage: 

a{r) = 1 \t\ 1+i3 , -oo < r < n, (6) 

where the inflation index /3 is an model parameter describing the expansion history during inflation. 
The special case of /3 = —2 corresponds the exact de Sitter expansion. However, both the model- 
predicted and the observed results indicate that the value of (3 could differ slightly from —2. 
The preheating stage : 

a(r) = a z \r - r p | 1+/3s , n < r < r„ (7) 

where the parameter (3 S describes the expansion behavior of the preheating stage from the end of 
inflation to the happening of reheating process followed by the radiation-dominant stage. In some 
literatures [41, 49], j3 a is set to be 1 , however, we take (3 S as a free parameter in the paper. 
The radiation-dominant stage : 

a(r) = a e (r - r e ), r s < r < r 2 . (8) 

The matter-dominant stage: 

a(r) = a m (r - r m ) 2 , r 2 < r < r E - (9) 

The accelerating stage up to the present time r [7]: 

o(r) = l H \r - r a |~ 7 , r E < r < r , (10) 

where 7 is a ^-dependent parameter, and is the energy density contrast. To be specific, we 
take 7 ~ 1.97 [50] for Q\ = 0.73 [34] in this paper. It is convenient to choose the normalization 
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l T o — T a\ — I; i- e -5 the present scale factor a(r ) = Ih- From the definition of the Hubble constant, one 
has Ih = ^/Hq, where H = 100 h km s _1 Mpc _1 is the present Hubble constant. We take h ~ 0.704 
[34] throughout this paper. Supposing (3 and (3 S are model parameters, all the constants included 
trough Eq.(6) to Eq. (10) can be fixed by the continuity of a(r) and a'(r) at the four given joining 
points Ti, t s , r 2 and t e , if one knows the increases of the scale factor of various stages, i.e., the 
definite values of Ci = a{r s )/a{n), ( s = a(r 2 )/a(r s ), C2 = a(T E )/a(r 2 ), and ( E = a(r )/a(T E ). 
The spectrum of RGWs h(k,r) is defined by 



2 , , N dk 



(^(r,x)%(r,x))=y o h 2 (k,r) T , (11) 

where the angle brackets mean ensemble average. The dimensionless spectrum h(k, r) relates to the 
mode /ifc(r) as [47] 

h(k,r) = ^k 3 / 2 \h k (r)\. (12) 

7T 

The one that we are of interest is the present RGWs spectrum h(k, tq). The characteristic comoving 
wave number at a certain joining time r x is give by [41] 

K = k(r x ) = YJW - y (13) 

After a long but simple calculation, it is easily to obtain kn = 27r7 and the following relations: 

kE _ ^ 2 _ A ks _ /- fc l _ f->A\ 

k H k E k 2 k s 

In the present universe, the physical frequency relates to a comoving wave number k as 

^ ~ 27ra(r ) ~~ 2irl H ' ( ^ 15 ' ) 

The present energy density contrast of RGWs defined by Qgw = (Pg) / Pc where p g = ^kG^^o ^ s 
the energy density of RGWs and p c = 3Hq/8ttG is the critical energy density, is given by [3, 5] 

rf upper 

V GW = / tog(f)j, (16) 

J flow J 

with 

n g (f) = ^hl(f)(pl 2 (17) 



H, 



being the dimensionless energy density spectrum. We have used a new notation, h c (f) = h(f, To) / \/2, 
called characteristic strain spectrum [5] or chirp amplitude [51]. The lower and upper limit of integra- 
tion in Eq.(16) can be taken to be f\ ow ~ f E and f upper — fi, respectively, since only the wavelength 
of the modes inside the horizon contribute to the total energy density. 



III. THE INCREASES OF THE SCALE FACTOR 



For the simple ACDM model, the late-time acceleration of the universe is well know. One easily 
has ( E = 1 + z E = (fiA/^m) 1 ^ 3 — 1-4, where z E is the redshift when the accelerating expansion 
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begins. The increase of the scale factor duration of the matter-dominated stage can also be obtained 
straightforwardly, £ 2 = ^^^^ = (1 + z eq)CE l with z eq = 3240 [34]. However, the histories of the 
radiation-dominated stage and the preheating stage are not known well. Recently, Mielczarek [40] 
proposed a method to evaluate the reheating temperature, Trh, under the frame of the slow-roll 
inflation model with a quadratic potential V(<t>) = \m 2 (f) 2 combing the observations from WMAP. 
Using this method, ( s and £1 can be determined subsequently with the evaluation of Trh [41]. In 
this paper, we consider a more general power-law form of the potential, V(<f)) = A0 n , where A is a 
constant. For this general form of V((j)), it is hard to obtain the analytic expression of the energy 
density of the universe at the end of inflation, and in turn, it is hard to obtain the temperature 
of reheating analytically. Hence, we will take some reasonable values of Trh constrained by CMB 
observations [43]. 

Firstly, we discuss the value of ( s . After reheating, the universe is filled with the relativistic plasma, 
which undergoes a adiabatic expansion as long as the entropy transfer between the radiation and 
other components can be neglected. The adiabatic approximation leads to the conservation of the 
entropy, i.e., dS = 0. It implies sa 3 =const, where the entropy density s of radiation is given by 

s = %9sT\ (18) 

Here, g s counts the effective number of relativistic species contributing to the radiation entropy. 
Another similar quantity g, counting the effective number of relativistic species contributing to the 
energy density of radiation, relates to energy density: 

(19) 

The behavior of g and g s with different energy scale were demonstrated in [46]. At the energy 
above ~ 0.1 MeV, one has g = g s . Moreover, at the energy scales above ~ 1 TeV, g = 106.75 
in the standard model, and g ~ 220 in the minimal extension of supersymmetric standard model, 
respectively. On the other hand, at the energy scales below ~ 0.1 MeV, g = 3.36 and g s = 3.91 
respectively. According to the conservation of entropy, one can easily gets the increase of the scale 
factor from the reheating till the recombination [40], 

dree TrH / 9*s 

a(r s ) T rec 

where a rec and T rec stand for the scale factor and the temperature at the recombination, respectively. 
g* s and g* a count the effective number of relativistic species contributing to the entropy during the 
reheating and that during recombination, respectively. As discussed in [43], the lower band of the 
reheating energy scale is 17.3 TeV constrained by the observed scalar power spectrum of CMB at 
95% of the confidence limit. Thus, in this paper we assume g* s ~ 200 eclectically, which was also 
employed in [43]. On the other hand, one has g+ s = 3.91 including the contributions of the effective 
number from photons and three species of massless neutrinos to the radiation entropy during the 
recombination, since the energy scale at the recombination T rec = Tcmb(1 + z rec ) ~ 10~ 7 MeV. 
Under the assumption of g* s ~ 200, the lower bound of Trh ~ 6 ■ 10 3 GeV was obtained [43]. On 
the other hand, gravitinos production gives an upper bound [52]. For instance, in the framework of 
the Constrained minimal supersymmetric standard model [44], the upper bound of Trh was found 
that Trh ^5 a few xlO 7 GeV from over-production of 6 Li from bound state effects, and moreover, 
Trh can be relaxed to < a few xlO 8 GeV when a more conservative bound on 6 Li/ 7 Li was used. 
However, if one does not consider the gravitinos production problem, the most upper bound of Trh 
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could be up to < 3 ■ 10 15 GeV coming from the energy scale at the end of inflation [43]. Based on 
Eq. (20), one easily obtain 

c= gw^. = ibh (g»y (21) 

a rec a{r s ) iCMB(i + z eq ) \g*sj 

where we have used T rec = T CM b(1 + z rec ). With T CM b = 2.725 K = 2.348 • 1(T 13 GeV [34], one has 
( s ~ 5 x 10 16 for example. Secondly, we discuss the evaluation of (±. First of all, we briefly recall 
the slow-roll inflation model. For slow-roll inflation, the evolution is described by the usual slow-roll 
parameters [53]: 



m? Pl fV'Y _m 2 Pl V" 



16tt V V J 8tt V 



(22) 



which are required to be much small than unity for the slow-roll approximation to be valid, e 
approaches to unity at the end of inflation. When the slow-roll conditions are satisfied, inflation 
continues keeping the Hubble rate nearly constant, and the primordial tensor power spectrum and 
the scalar power spectrum are respectively given as [48, 51]: 

A'for^-f^Y, (23) 

A^(fc,T,)«-f— V, (24) 
ire \mpij 

where if* is the Hubble rate during inflation, and r* stands for the moment when the /c-mode exits the 
horizon. On the other hand, based on the observations of CMB, the present scalar power spectrum 
can be expanded in power laws, 

A* (A:) = Al(k ) if]" , (25) 



k 



k x ns " 1 



Ai(k) = Ai(k ) ( - 1 , (26) 

where A^(fc ) and A^(k ) are the power spectrum of the tensor perturbations and curvature per- 
turbations evaluated at the pivot wave number k% = fco/ a ( r o) = 0.002 Mpc -1 [34], respectively. 
Furthermore, under the slow-roll approximation, at the pivot wave number fc the spectral parame- 
ters are given by [53] 

rat ^ -2e, (27) 
n s ~ 1 - 6e + 2r), (28) 

In general, the spectral indices n t and n s are k— depedent, described by the running parameters 
at = dn t /d\nk and a s = dn s /d\nk, respectively [47, 53-55]. However, at and a s are only second 
order small quantities. Moreover, if one uses the quantum normalization (see below) as the initial 
condition for the generation of RGWs, a t should be exactly zero. On the other hand, as will be 
seen below, non-zero a s would induce an n s greater than 1, which make us difficult to evaluate 
the increase of the scale factor from the ko mode exiting the horizon during inflation to the end of 
inflation. Hence, in this paper we will simply set at = a s = 0. Note that Even though the value of nt 
is quite uncertain, n s can be well constrained by CMB [34] or BAO [56] . The ratio of the primordial 
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tensor power spectrum to the scalar power spectrum is defined as [48, 51] 



A^(fc, T*) 

based on Eqs.(23) and (24). Therefore, at the pivot number k , one has 

_ Al(ko,n) ^ Al(k ) 



16e, (29) 



(30) 



where Tj is the fco-mode exit the horizon during inflation. The approximation of the second equation 
in Eq.(30) accounts for that the pivot ko wave mode reentered the horizon a little earlier than the 
present time, and then has suffered a decay. Therefore, the ratio A^(fe )/A^(fco) can not exactly 
reflect the true value of r given by its definition, however, the deviation would be expected to be less 
than ~ 0.8% [41]. Hence, we will use this approximation when confront with the CMB observations. 
Furthermore, under this approximation, one has a simple relation: 

n t = 2/? + 4, (31) 

since the primordial spectrum of RGWs has a power-law form A(fco) — A(/co?T"j) oc k 2/3+4 [41]. 
WMAP 7 Mean [34] fixed A^(fc ) = A, = (2.43 ± 0.11) • 10~ 9 . Thus, the non-zero value of r implies 
the existence of gravitational wave background, which induced uniquely the B-mode polarization 
of CMB [57]. At present only observational constraints on r have been given [33, 34]. The upper 
bounds of r are recently constrained [34] as r < 0.24 by WMAP+BAO+if and r < 0.36 by WMAP 
7 only for a s = 0, and r < 0.49 for a s ^ by both the combination of WMAP+BAO+ifo an d 
the WMAP 7 only, respectively. Furthermore, using a discrete, model-independent measure of the 
degree of fine-tuning required, if 0.95 < n s < 0.98, in accord with current measurements, the tensor- 
to-ratio satisfies r > 10~ 2 [58]. Therefore, one can normalize the RGWs at k = ko using Eq. (30), if 
r can be determined definitely. 

As analyzed by Mielczarek [40], for the pivot wave number fcg, the total increase of the scale factor 
from the mode exit the horizon during inflation up to the present time can be evaluated as 

-£/* 
kn 



Ctot - "TjT- (32) 



•o 



Due to Eqs. (24) and (26), one has 



ire \ mpi 



A^(fco), (33) 



where the approximation A^(A; ) ~ A^(A; ,Ti) was used. Taking the form V(<f>) = X(f> n , one can 
easily have a relation: 

nil — n s ) , . 

e = -A y~ 34 

2(n + 2) v ' 

from Eqs. (22) and (28). Plugging Eqs. (33) and (34) into Eq. (32) gives 



^ a f&-^w (35) 

On the other hand, if we assume the universe did a quasi-de Sitter expansion(/3 w —2), the increase 
of a scalar factor from the moment of ko mode exiting the horizon during inflation to the end of 
inflation is give by 

d = e N , (36) 
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where N is the e-folding number, which can be estimated as 



N 



Concretely, for V((j)) = A0 n , one can get 



v{4>) 



N 



n + 2 

2(1 - n y 



(37) 



(31 



with the help of Eqs. (22) and (28). So, if n = 2, Eq.(38) 
Plugging Eq. (38) into Eq. (36), and using the identity 

Ctot = OC1CC2CB; 

one can easily obtain the complete expression of (1. 



reduces to the result shown in [40]. 



(39) 



Ci 



mpi 



7rA^(fc )(l - n s 
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2 n + 2 



1/2 



T, 



CMB 



RH 



9*8 

9*s 



1/3 



— exp 



n + 2 



2(1-71, 



(40) 



One can examine that, for tj = 2, the above expression reduces to Eq. (11) in Ref.[41] after using 
Eq.(7) in the same reference. In the following, let us see the reasonable range of the index n 
constrained by both theories and observations. As well known, at the end of inflation, the scalar 
field <p oscillates quickly around some point where V(<f)) has a minimum. In the limit that the 
oscillation rate is much greater than Hubble expansion rate H, and ignoring the coupling between 
the scalar field <j) and other components, it is found that [42] the scalar field oscillations behave like 
a fluid with p = wp, where the average equation of state w depends on the form of the potential 
V ((/)). For V{<j>) = \4> n , one has 



w 



n-2 
n + 2 



(41) 



and p decreases as a ~ 6n /( n + 2 ). i n particular, n = 2, one has w = and p oc a~ 3 , which imply 
a matter-dominant like expansion of the preheating stage [49]. Adding the consideration of the 
coupling between the scalar field and the resulting relativistic particle creation, Martin and Ringeval 
[43] verified the relation (41) using a numerical method, and it was found that the average w never 
deviates from zero exceeding 8%. From theoretical consideration, one should have w < 1 to satisfy 
the positivity energy conditions; while w > —1/3 to make sure the inflation must stop and the 
preheating stage begins. Due to Eq. (41), the condition — 1/3 < w < 1 leads to n > 1. On the 
other hand, Martin and Ringeval [43] firstly gave a constraint on n based on the CMB observations, 
n < 2.1. Therefore, based on both the theories and observations, the index n is constrained to be 



1 < n < 2.1. 



(42) 



Note that, there is a relation between n and f3 s . According to the energy conservation equation and 
the Friedmann Equation, 



a 

p + 3-p(l + w) 

'a\ 2 8ttG 
a I 3 



0. 



(43) 
(44) 
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FIG. 1: The relation between r and n for the fixe value of n s = 0.96, n a — 0.97 and n a = 0.98, respectively. 

one can easily obtain a oc £ 2 /(3+3w) ^ r 2/(i+3w)_ TJsing Eqs. (7) and (41), and allowing for p oc 
a _ 6n/(n+2 \ one has 

A — n 

ft = 25^- (48) 

Then, in principle, the expression of (j in Eq. (40) can be rewritten as a function of /3 S . From the 
combination of Eqs. (42) and (45), one finds that, n > 1 leads to f3 s > —0.5 and n < 2.1 leads to 
/3 S > 0.86, respectively. Based on the range of n (or (3 S ) discussed above, we try to constrain some 
parameters combining with CMB observations. 

IV. PARAMETERS CONSTRAINTS FROM OBSERVATIONS 

As shown in the previous section, many parameters are dependent on the value of n s . Seven-year 
WMAP Mean [34] gives n s = 0.967 ±0.014, and n s = 0.982±g;g?g when one also considers the tensor 
mode contributions to the anisotropies of CMB. Moreover, the combination WMAP+BAO+i^o Mean 
gives n s = 0.968 ± 0.012, and n s = 0.973 ± 0.014 when the tensor mode contributions are included. 
Independently, SDSS III predicts n s = 0.96 ± 0.009 [56]. As can be seen in Eq. (40), ( is sensitively 
dependent on n s , and in turn one can expect that the spectrum of RGWs also depend sensitively 
on n s in the very high frequencies. Therefore, for a general demonstration, we consider the cases: 
n s = 0.96, 0.97 and 0.98, respectively. 

Firstly, let us constrain the tensor-to-scalar ratio r. According to Eqs. (22) and (29), it is straightly 
to get 

l-n s ). (46) 



n + 2 

We show this relation in Fig.l. One can see that r increases slowly with n. r lies in (0.11, 0.16), 
(0.08, 0.12), and (0.05, 0.08) for n s = 0.96, n s = 0.97, and n s = 0.98, respectively. Similarly, from 
Eqs. (27) and (31), one has 

71 

/3 = -2- — - r (l-n s ), (47) 
2(n + 2) v ; ' v ; 

which is shown in Fig. 2. The parameter (3 is constrained in the range of (—2.007,-2.010), 
(-2.005,-2.008), and (-2.003,-2.005) for n s = 0.96, n s = 0.97, and n s = 0.98, respectively. 



10 




FIG. 2: The relation between /3 and n for the fixe value of n 3 = 0.96, n 3 = 0.97 and n s — 0.98, respectively. 




Therefore, the range of n in Eq. (42) leads to very tight constraints on r and (3, which are limited 
in very narrow ranges with definite value of n s . 

Now, let us see the increase of the scale factor during preheating stage d, which is expressed in 
Eq.(40). We plot it in Fig. 3 as a function of n with definite values of Trh. Allowing for the expansion 
of the universe, one would expect that Ci > 1- As can be seen in Fig.3, the cases of n s = 0.96 and 
0.97 can make sure well the resultant £i being much larger than 1, however, the case of n s = 0.98 
can not be compatible with the fact £i > 1 in the whole range of n shown in Eq.(42). If n s is 
determined well to be as high as 0.98, it will give very tight constraints on n. Concretely, n < 1.7 
and n < 1.3 for Trh = 10 4 GeV and Trh = 10 8 GeV, respectively. What we are more interesting 
are the characteristic frequencies given by Eq.(15). With the help of Eq. (14), one can easily get the 
characteristic frequencies: f H = H ~ 2.28 • 10~ 18 Hz, f H = H ~ 2.28 ■ 10~ 18 Hz, f E ~ 1.93 • 10~ 18 
Hz, / 2 ~ 9.3- 10~ 17 Hz. The value of f s depends linearly on Trh- For instance, f s ~ 4.54- 10 -3 Hz for 

T RH = 10 4 GeV and f s ~ 45.4 Hz for T RH = 10 8 GeV, respectively. Since f x = fsCl^ 7 , it depends 
on the values of n (or (3 S ), n s and Trh- It is worth to be study in detail, since the characteristic 
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FIG. 4: The upper limit frequency /i as a function of n for different combinations of n s and Trh. The dotted parts for n a = 0.98 
are constrained by the condition C,i > 1 shown in Fig. 3. 



frequency f\ is approximately the highest frequency of RGWs. The modes whose frequency higher 
than fi decay with the expansion of the universe [2, 3]. We plot f\ as a function of n with definite n s 
and Trh in Fig.4. One can see that the behaviors of f\ along with the increasing n are quite different 
for different values of n s . For the case of n s = 0.98, we plotted j\ using dotted lines for the values 
of n constrained by (i > 1 which are shown in Fig.3. In the part of large n, f\ is larger for smaller 
values of n s . On the other hand, in the limit of n — > 1, fi becomes a fixed value independent on n s , 
and moreover, the asymptotic fixed f\ is larger for a larger value of Trh- It is easy to understand if 
one has found that /3 S — > +oo as n — > 1 from Eq.(45) which leads to /i — > f s . The value of f\ should 
be below the constraint from the rate of the primordial nucleosynthesis, /^3x 10 10 Hz [2]. When 
the acceleration epoch is considered, the constraint becomes f\ ~ 4 x 10 10 Hz. This will in turn give 
some constraints on n, n s and Trh- 

As analyzed in our previous work [41], when the quantum normalization for the generation of 
RGWs during inflation is employed, one has 

A n (k y/ 2 = SV^hiHoC^ C^cfcT (j±Y, (48) 

where Ipi = \fG is the Planck length. In Eq.(48), there are totally six parameters: r, f3, /3 S , n, 
Trh and n s . However, among them only three are independent, due to Eqs.(45), (46) and (47). We 
show the Trh — relation with definite values of n s in Fig. 5. First of all, we define the range of 
6 • 10 3 — 10 8 GeV as Region I; while the range of 6 • 10 3 — 3 • 10 15 GeV as Region II. It is found 
that, under the condition of quantum normalization, n s = 0.96 and n s = 0.98 can be ruled out, 
since the resultant Trh outside Region II. If one consider the gravitinos production problem, the 
case n s = 0.97 would also be ruled out, since the resultant Trh outside Region I. However, the 
resultant Trh given by n s = 0.966 lies well inside Region I for the whole range of /3 given by Eq. 
(47). Moreover, as shown in Fig. 5, the quantum normalization will give a little tighter constraints 
on the range of (3 for n s = 0.967 and n s = 0.968. Note that, these results are based on the validity 
of quantum normalization, however, it is not the unique initial condition. Let us make a comparison 
with the previous results in [41]. Taking n s = 0.966 for example, n = 2 leads to Trh — 3.4 ■ 10 6 
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FIG. 5: The relation between Trh and /3 based on Eq.(48) due to the condition of quantum normalization. 



Without quantum normalization 




f[Hz] 

FIG. 6: The energy density spectra of RGWs for the fixed p s = 1 with different combinations of n s and Trh, without considering 
the quantum normalization. 

GeV; while Trh — 2.8 ■ 10 12 GeV shown in Fig.l in Ref.[41]. Hence, the discrepancy of Trh, at six 
orders of magnitude, indicates that the quantum normalization may be not a good initial condition. 
However, one should also keep in mind that we have used many approximations, which would also 
contribute a lot to the discrepancy of Trh discussed above. Note that, if one does not consider 
quantum normalization, the zero point energy should be removed or else the cosmological constant 
would be 120 orders of magnitude larger than observed. Some effective methods [59] have been 
pointed out. In next section, we will demonstrate the spectra of RGWs with and without quantum 
normalization respectively. 
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Without quantum normalization 




f[Hz] 

FIG. 7: The energy density spectra of RGWs for the fixed Trh = 10 6 GeV with different combinations of n a and f3 s , without 
considering the quantum normalization. 

V. THE SPECTRA OF RGWS AND THEIR DETECTION 

In this section, we demonstrate the energy density spectra of RGWs with reasonable values of the 
parameters and discuss the detection due to the current running and planned gravitational wave 
detectors. 

As discussed in the previous sections, there are many parameters involved in the spectrum of 
RGWs. They are n s , n, r, (3, j3 s , and Tr H . However, among them only three are independent due to 
Eqs. (45)-(47). Furthermore, n s has been constrained well from observations of CMB, BAO and H . 
Since the spectrum of RGWs in the high frequencies extreme sensitively depends on n s , we discuss 
two cases of n s = 0.968 and n s = 0.973, respectively, based on the combination of WAMP+BAO+ifn 
[34]. In the following, we regard (3 S and Trh as parameters, and choose some representative values 
of them since they have large uncertainties. In order to give a complete discussion, we will consider 
the spectra of RGWS both with and without quantum normalization. As analyzed in Sec. Ill, Trh 
is constrained to be Trh ~ 10 4 — 10 8 GeV, and /3 S is limited to be larger than 0.86. Firstly, let us 
see the case of no quantum normalization. Setting f3 s = 1, Fig. 6 shows the energy density spectra 
of RGWs, Ogw(/), with different values of n s and Trh- One can see that, all the Ogw(/) nearly 
overlap each other in the low-frequencies. This is because the spectrum for / < f s is only related to 
r and (3 [41], and, moreover, the differences of r and (3 are very small between the case n s = 0.968 
and the case n s = 0.973 due to Eqs. (46) and (47). Explicitly, one has r = 0.128, (3 = —2.008 for 
n s = 0.968 and r = 0.108, (3 = —2.007 for n s = 0.973, respectively. However, in the part of high 
frequencies, Oqw(/) exhibits different properties for different combinations of n s and Trh- On one 
hand, for the same value of Trh, the spectrum ficw(/) with n s = 0.968 and that with n s = 0.973 
have the same "turning point" from which Ogw(/) decreases rapidly with the increasing frequency, 
and the "turning point" is just f s which is only dependent on Trh- Moreover, the decreasing slope 
of the logarithm of the two spectra for / > f s are nearly the same since f^Gw(/) ^ f 4+2l3 ~ 2/3a [41] 
which is reduced to Qqw(/) ^ f~ 2 f° r ~ — 2 and f3 s = 1. However, the Ogw(/) with a smaller 
n s has a larger upper limit frequency fi which responds to a lower amplitude of Ogw(/)- On the 
other hand, for the same value of n s , the G w(/) with a higher T RH leads to not only a larger f s 
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With quantum normalization 




f[Hz] 

FIG. 8: The energy density spectra of RGWs under the condition of the quantum normalization. 

1 /2 

but also a larger fi since fi oc T R y for (3 S = 1 which can seen from the combination of Eqs. (14), 
(21) and (40). Fig.7 shows the energy density spectra of RGWs for the fixed value T R h = 10 6 GeV. 
One can see that a larger (3 S leads to a steeper slope of the logarithm of f2cw(/) and a smaller fi 
for the same values of n s . In a word, Trh determines the value of f s , (3 S determines the slope of the 
logarithm of ^gw(/) f° r the fact that (3 ~ —2, and fi depends on all the three parameters especially 
n s . Secondly, let us consider the case of the quantum normalization. Due to the resultant Eq. (48), 
among the three parameters Trh, /3 s and n s only two of them are independent. Taking n s and Trh 
as parameters, Oqw(/) with some combinations of the two parameters are plotted in Fig.8. 

Below, let us discuss the detection of RGWs using the ongoing and planned gravitational detectors 
which are sensitive at different frequency bands. As shown in Fig.6-Fig.8, the differences of the 
spectra of RGWs with different parameters are only significant in high frequency parts. Hence, we 
just take a characteristic combination of the parameters n s = 0.968, f3 s = 1 and Trh = 10 6 GeV for 
demonstration. As a conservative evaluation, in Fig. 9 we show the strain amplitudes, h c (f)/yfj [5] 
of RGWs confronting the strain sensitivity curves of various gravitational wave detectors including 
the complete PPTA [25] and SKA [60] using the pulsar timing technique, and the space-based laser 
interferometers such as eLISA [19], BBO [20, 21], and the Fabry- Perot DECIGO [22]. One can see 
that, RGWs under the frame of the slow-roll inflation with a potential V(<f)) oc <p n are quite promising 
to be detected by the future SKA, eLISA, BBO and DECIGO. As seen from Fig.6-Fig.8, fl GW (f) with 
different parameters have different properties in high frequencies. It would be interesting to discuss 
the detection of RGWs in high frequencies. In Fig. 10, we plot the characteristic amplitude of RGWs 
with different parameters and conditions compared to the instrumental noise, a/ fS n (f), of BBO, the 
ultimate DECIGO [23], the second generation ground-based laser interferometers AdvLIGO [11], and 
the third generation ET [17]. The parameters and conditions of RGWs are listed in Table I. S n (f) 
is the normal one-side noise spectrum of detectors. As can be seen in Fig. 10, even though AdvLIGO 
and ET are hard to catch the signals of RGWs, BBO has the potential to distinguish RGWs with 
different parameters or different conditions in the frequency band 10~ 2 — 10° Hz. Furthermore, 
the ultimate DECIGO has the capability to distinguish them more easily. Thus, the future BBO 
and DECIGO detections provide an important tool not only determining the parameters but also 
examining the validity of the quantum normalization when RGWs were generated during inflation. It 
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FIG. 9: The strain of RGWs with different parameters for n B = 0.968, f3 s = I and Trh = 10 6 GeV confronts against the current 
and planed GW detectors. The sensitivity curves of PPTA and SKA using pulsar timing technique are taken from Refs.[25] 
and [60], respectively. The curve of BBO is generated using the online "Sensitivity curve generator" [19] with the parameters 
in Table II of Ref.[20] and Table I of Ref.[21]. The curve of DECIGO is taken from Ref.[22]. 
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FIG. 10: The characteristic amplitude of RGWs in high frequencies confronting against the instrument noise y/ fS n (f) of BBO, 
Ultimate DECIGO [23], AdvLIGO and ET. 

is worth to point out that, at frequencies lower than 10 -2 Hz the signals of RGWs are contaminated 
by the confusion noise produced by galactic binaries [61]. Hence, we should focus on the frequencies 
higher than 10~ 2 in order to distinguish various spectra of RGWs. 

VI. CONCLUSIONS AND DISCUSSIONS 

In the frame of the slow-roll inflation with a power-law form V = X(p n , we calculated the analytic 
solutions of RGWs. In the narrow range 1 < n < 2.1, the tensor-to-scalar ratio r and the inflation 



16 



TABLE I: The definitions of h c with different parameters. "N" stands for "No" meaning that the condition of the quantum 
normalization is not considered; while "Y" stands for "Yes" meaning that the condition of the quantum normalization is 
considered. 



h c 


n s 


ft 


Trh 


quantum normalization 


hd 


0.968 


1 


10 4 GeV 


N 


h c 2 


0.968 


1 


10 6 GeV 


N 


h C 3 


0.968 


4 


10 6 GeV 


N 


h c 4 


0.973 


1 


10 6 GeV 


X 


h C 5 


0.968 


1 


2 ■ 10 11 GeV 


Y 



expansion index (3 are both tight limited to lie in narrow ranges for a given value of the scalar spectral 
index n s . Concretely, r lies in (0.11, 0.16), (0.08, 0.12), and (0.05, 0.08) for n s = 0.96, n s = 0.97, 
and n s = 0.98, respectively; while [5 lies in the range of (-2.007,-2.010), (-2.005,-2.008), and 
(—2.003, —2.005) for n s = 0.96, n s = 0.97, and n s = 0.98, respectively. Moreover, the preheating 
expansion index j3 s is constrained to be /3 S > 0.86. We found that the spectrum of RGWs in high 
frequencies depends on the parameters n s , f3 s and Trh- Explicitly, Trh determines f s where the 
flat RGWs spectrum decreases suddenly. (3 S determines the decreasing slope of the logarithm of the 
spectrum. Whereas, the upper limit frequency f\ is dependent on all the three parameters n s , (3 S 
and Trh- Besides, the quantum normalization for the generation of RGWs also affect the spectrum 
of RGWs in high frequencies. 

Among the current and planed GW detectors, SKA using the pulsar timing technique and the 
space-based interferometers eLISA, BBO and DECIGO are promising to catch the signals of RGWs. 
Furthermore, BBO and DECIGO have the potential not only to distinguish the spectra with different 
parameters but also to examine the validity of the quantum normalization. Therefore, RGWs could 
become the most important tool to know the physics occurred in the very early Universe such as 
the inflation and reheating process. Even though we chose a series power-law form potential of the 
scalar field, as shown in [42], a polynomial form of the potential will be dominated by the lowest 
power of in V. In this case, the conclusion is not substantially modified. In our previous work [41], 
we got the r — (3 relation for a particular potential V = ^m 2 (f) 2 . However, for the more general case 
V = \<p n , it is hard to obtain a complete analytic solution of Trh and in turn the increases of the 
scale factor ( s and (j. Therefore, in this paper, we set a series reasonable values of Trh as additional 
parameters. The determination of n s is very important, since it is sensitively affect our results. The 
future CMB experiments such as the Plank satellite [35], the ground-based ACTPol [36] and the 
planned CMBPol [37] will help us to determine the more convincible value of n s . Therefore, one can 
expect accordingly that the spectrum of RGWs would be known better. 

In principle, our analysis is valid for the slow- roll inflation with other forms of the potential V(<j)). 
However, for some particular forms of V(<f)), it would be difficult to get the analytic result of d 
as a function of the parameters included in V{<p). Moreover, one could not effectively constrain 
the parameters in V(<fr). However, one can still calculate £i numerically according to the whole 
calculating processes presented in Ref.[40], and then calculate the spectra of RGWs accordingly. 
More general inflationary models other than the slow-roll inflation and the slow-roll inflation with 
other forms of V((f)) would be studied in our future work. 
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